Sergio Oller, Maritzabel Molina, Juan Carlos Vielma, Xavier Martínez, Alex H. Barbat 
y Liz Nallim. Capítulo 11 -Análisis numérico de la reparación y refuerzo de estructuras 
con FRP. Omnia Publisher, 2014. p. 261-291. Aplicaciones avanzadas de los 


materiales compuestos en la obra civil y la edificación ISBN: 978-84-941872-9-2, 
Diposito legal: B-12444-2014, DOI: 10.3926/0ms.209. 


11 


Análisis numérico de la reparación 
y refuerzo de estructuras con FRP 


Sergio Oller, Maritzabel Molina, Juan Carlos Vielma, 
Xavier Martínez, Alex H. Barbat y Liz Nallim 


Resumen En este capítulo se presenta un procedimiento para evaluar la magnitud del daño 
local y el daño global en estructuras sometidas a acciones estáticas y dinámicas, con especial 
énfasis en el problema sísmico. Además de la formulación para la evaluación del daño, se 
introduce el concepto de refuerzo y reparación estructural mediante laminados compuestos 
de matriz epóxica con refuerzo de fibras de carbono. Para esta finalidad, se menciona la 
utilización de la teoría de mezclas para componer un material compuesto a partir de sus 
componentes básicas. También se evalúa el daño en estas estructuras reforzadas y/o reparadas 
y se comenta la influencia de estas mejoras en la valoración del daño global de la estructura. 


11.1. Introducción al comportamiento no lineal de las es- 
tructuras y valoración del estado de daño 


En este capítulo se presenta una descripción metodológica que permite evaluar el daño 
global en estructuras de hormigón armado sometidas a acciones estáticas/dinámicas. Prime- 
ramente se define el estado de daño local de un punto de la estructura, al estado mecánico que 
resulta luego que este inicie un proceso no lineal disipativo provocado por diversos mecanis- 
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mos internos (pérdida de resistencia por efecto de la degradación de rigidez, por incremento 
de deformaciones irreversibles, por fatiga cíclica, etc.). A partir del concepto de daño local 
en un punto del sólido, se deduce por extensión el índice de daño global estructural [6]. 


En esta descripción metodológica que aquí se presenta, se utiliza la técnica de elementos 
finitos para aproximar el campo de desplazamientos y a partir de él obtener el estado de de- 
formación con los que se calculará el correspondiente estado de tensión en cada punto de la 
estructura. La formulación es general y se puede utilizar en estructuras tridimensionales, bi- 
dimensionales y elementos estructurales de láminas y barras. En los elementos estructurales 
se utiliza una formulación apropiada para placas de Reissner-Mindlin [37] y barras Timos- 
henko [37] subdivididos en capas entrecruzadas (fibras en las vigas). Esto permite obtener los 
estados tensionales elástico predictivo en cada punto de la estructura y a partir de ellos utili- 
zar un modelo constitutivo formulado en tensión-deformación para obtener el correcto estado 
tensional del punto. Esta forma de trabajar permite incluir cualquier fenómeno mecánico que 
acerque más la simulación numérica a la realidad del problema (plasticidad, degradación de 
rigidez, problemas viscosos, etc.). 


En otros trabajos ya publicados por los autores [1, 2, 7, 27], se ha desarrollado una metodo- 
logía de simulación del daño ocurrido en estructuras de hormigón armado durante terremotos. 
Dicha metodología se apoyaba en un modelo constitutivo de degradación isótropo [22], que 
derivaba de la teoría de Kachanov [8]. En aquellos trabajos, se entiende por daño el grado de 
degradación estructural -pérdida de rigidez- y se caracteriza por un índice escalar normaliza- 
do que varía entre O < d < 1, según el material esté en estado virgen o totalmente degradado 
sin resistencia residual, respectivamente. 


En este trabajo, se extiende dicho concepto con la finalidad de que este índice pueda incluir 
otros fenómenos disipativos como la plasticidad y la degradación a la vez, permitiendo un 
mejor acercamiento al concepto de deterioro en un punto de la estructura. Este índice más 
general, basado en la normalización de fuerza residual no equilibrada, preserva la definición 
del índice de daño local con la misma objetividad que se ha propuesto en trabajos previos 
[1, 2, 7, 27], es decir, como un escalar acotado entre cero y uno que permite definir el estado 
del punto desde no dañado a totalmente dañado, respectivamente. 


11.2. Introducción a la simulación numérica de estructuras 
de hormigón armado. Enfoque conceptual 


Las estructuras de hormigón simple y también las reforzadas con acero (hormigón armado) 
y con FRP (Fiber Reinforced Polymer), están conformadas por un material compuesto de más 
de un material simple (matriz: cemento más agregado pétreo, refuerzo: acero y/o FRP). Estas 
estructuras tienen un comportamiento inicial elástico hasta que algún material componente 
del compuesto supere el correspondiente umbral de resistencia, instante en que se inicia el 
comportamiento no lineal de dicho componente, induciendo también al compuesto a un com- 
portamiento no lineal. De esta manera, si algún punto de la matriz de hormigón supera el 
umbral de resistencia definido en mecánica por alguna función de discontinuidad (por ejem- 
plo Mohr-Coulomb, Drucker-Prager, etc.), entra en un estado de micro fractura y posterior 
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coalescencia que luego conduce a macro fracturas. Simultáneamente el acero, por condicio- 
nes de compatibilidad limitada (estado de comportamiento conjunto entre los dos materiales), 
continúa con un proceso de comportamiento elástico con posterior deslizamiento incompati- 
ble respecto de la matriz de hormigón. De mantenerse la solicitación monótona creciente, el 
acero podría superar su umbral de discontinuidad (von Mises) y alcanzar un comportamien- 
to plástico. Esta diversidad de comportamiento entre los materiales componentes (matriz de 
hormigón y refuerzo de acero y/o FRP), que coexisten inicialmente en un estado elástico con 
movimientos compatibles, pierden luego dicha compatibilidad y se produce un deslizamiento 
relativo entre los materiales componentes adyacentes. Esta diversidad de comportamiento no 
lineal en cada material y en su interfaz resulta muy difícil de representar mediante una única 
ley constitutiva. En tal sentido, en este capítulo se muestra un procedimiento general basado 
en la mecánica de medios continuos y en técnicas numéricas, que permiten abordar el com- 
plejo comportamiento de los materiales compuestos en general y en particular el hormigón 
armado. 


Las bases conceptuales de este trabajo se asientan en la modelización constitutiva de cada 
simple material componente, que luego participará en el tratamiento de composición del ma- 
terial compuesto mediante la teoría de mezclas de sustancias básicas [26]. No obstante, estas 
dos formulaciones mecánicas permiten abordar gran parte de este complejo problema, sien- 
do también necesario completar la formulación incluyendo la anisotropía de los materiales 
componentes. Este es el caso del acero, cuyo comportamiento es isótropo, pero su utilización 
en forma de barras de refuerzo hace que su comportamiento fenomenológico resultante sea 
ortótropo uniaxial en la dirección que está orientada la barra. 


Con base en lo antes mencionado, es importante notar que gran parte del acierto en la res- 
puesta del material compuesto, comienza por una buena representación del comportamiento 
básico de cada material componente a través de adecuados modelos constitutivos. Son estos 
modelos no lineales, junto con el correspondiente a la interfaz entre los materiales componen- 
tes, los que inducen los cambios de comportamiento en el compuesto, pasando de la linealidad 
a la no linealidad. 


La no linealidad se manifiesta a través de una falta de equilibrio en la estructura global (Af, 
fuerza residual no equilibrada), como consecuencia de que algunos materiales componentes 
que integran la estructura, no pueden soportar los estados tensionales a los que se encuentran 
sometidos. Esto se puede formalizar en la siguiente ecuación de equilibrio dinámico en un 
cierto instante del proceso mecánico: 


Af = £”% (Ù, 1) 4£ (0,U, t) -£% (1) =0 (11.1) 


donde £”"** (Ü, t) = MÜ (t) es la fuerza másica o inercial que depende de la masa y la ace- 
leración Ü de cada punto. £(U, U, t) = DÙ (t) + KU (t) = fy O VÍN (x,y,z) dV es la fuerza 
interior que considera la rigidez K y la amortiguación estructural D, que actúan respectiva- 
mente sobre el desplazamiento U y la velocidad Ù de cada punto de la estructura. O = Y .k. O 
es la tensión desarrollada en un punto del sólido para cada componente c*%74 y N (x,y,z) es 
la función de forma o de interpolación local (MEP). Por último f®" (t) es la fuerza exterior o 
acción solicitante que actúa en el mismo tiempo 1 que las fuerzas antes mencionada. 
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El desequilibrio que se manifiesta en la ecuación 11.1, se elimina mediante una redistribu- 
ción del comportamiento de los puntos vecinos al que ha incursionado en el campo no lineal. 
Esta redistribución deberá conseguirse bajo un estado de deformación estructural compatible 
que desarrolle unas nuevas fuerzas interiores capaces de reequilibrar la estructura. Este su- 
puesto de convergencia hacia el equilibrio no siempre tiene solución, pues es posible que no 
se pueda alcanzar un estado de equilibrio compatible que equilibre la acción impuesta con las 
fuerzas interiores (véase ecuación 11.1 y Figura 11.3). 


Como puede verse en la representación esquemática de la Figura 11.1, en el grado de 
libertad correspondiente al empotramiento izquierdo de la viga biempotrada, el momento 
solicitante M®* no puede ser equilibrado (M®* — M'" = AM 4 0) por el momento interior 
M'" = f, 0zdA desarrollado en la sección transversal A por las tensiones interiores O. Esto 
suele ocurrir porque la resistencia del material no es ilimitada y las tensiones desarrolladas 
superan el umbral de discontinuidad elástica (f(0) — fy > 0), conduciendo a la aparición de 
micro fisuras y la posterior relajación de la tensión hasta cumplir con la capacidad límite del 
material f(0) — fy = 0. En el proceso de relajación y redistribución tensional, la estructura 
cambia su configuración cinemática, pasando de un giro Q; a otro @; + AQ; hasta alcanzar una 
nueva condición de equilibrio AM. Esta búsqueda de la solución compatible y equilibrada se 
consigue mediante el procedimiento de Newton-Raphson. Es decir, admitido el desequilibrio 
AM (t), esto es: 

AM(0) = M” (t) 44M" (0)-M“(0) 40 (11.2) 


ql 


Estado: 1 


M, 70:44, ,7+0 


0-1 


q2>q1 
Estado: 2 


M, 
M, #0< M, ; Mb, ,7+0; 6, >0, 


q3>q2 
Estado: 3 


M, 
M,=0<M,:; Ab, #0; 4,>0, 


Figura 11.1: Representación esquemática del procedimiento de Newton-Raphson para tratar el comportamiento no 
lineal de un grado de libertad rotacional de una estructura biempotrada. 


Se intenta a continuación eliminar el residuo AM (t) en el tiempo £, mediante la linealiza- 
ción en serie de Taylor de esta ecuación, truncada en la primera variación, 

ð [AM;(+)] 
dH 
— 
Ki(t) 


AM;-1(t) + A0¿=U > AQ; = —K;(t)AM;_1(£) (11.3) 


de donde resulta el incremento de giro AQ;, habitualmente denominado giro en la rótula plásti- 
ca, aunque no necesariamente tenga alguna relación con el conocido fenómeno de plasticidad. 


Ecuación de equilibrio dinámico y linealización del problema no lineal l 11.3 
En esta ecuación, K;(t) es la rigidez en el tiempo t, iteración de equilibro pa correspon- 
diente al grado de libertad rotacional del apoyo izquierdo y AM;_¡ (t) es la fuerza residual en 
el paso anterior. 


En la resolución general (véase sección 11.3), este procedimiento se aplica a todos los 
grados de libertad a la vez, resultando la ecuación de equilibrio no lineal mostrada en la ecua- 
ción 11.1. La solución de la misma mediante linealizaciones sucesivas (Newton-Raphson) 
[23, 37] se muestra en la Figura 11.2. 


11.3. Ecuación de equilibrio dinámico y linealización del 
problema no lineal 


La ecuación de equilibrio dinámico de un sólido discreto sometido a acciones externas 
variables en el tiempo, puede obtenerse directamente a partir de la Primera Ley de la Termo- 
dinámica [9, 10, 13]. A partir de los conocimientos previos sobre el método de los elementos 
finitos que se considera tiene el lector [37], puede escribirse la ley de conservación en la 
siguiente forma [23, 27], 

Fa 5 A A du; 

Joy: 1jdV = prindS+ | pbiiav — | pi 52 av (11.4) 

V S V V 
donde t; es la fuerza de superficie aplicada sobre el contorno S, (siendo t; = O;;n j, tal que O;; 
es el tensor de tensiones de Cauchy y n; el vector normal a la superficie S que envuelve el 
sólido); b; las fuerzas de volumen por unidad de masa; p = dM/0V la densidad de masa, 
siendo M la masa y V el volumen; v; = du¡/dt = ú; es el campo de velocidades (si t = cte. 
entonces la velocidad v; = ú; se transforma en un incremento temporal del campo de los 
desplazamientos, y la potencia introducida se transforma en el incremento temporal de tra- 
bajo introducido). La velocidad de deformación, ahora incremento temporal de deformación, 


puede escribirse como D;¡= {Li j} s= {Vù j} s= (E aN T que sustituida en la anterior, 
resulta el equilibrio de potencias en un sólido continuo. Además, basado en la técnica de 
elementos finitos, se aproxima mediante una función polinómica normalizada a la unidad 
N g(x, y,z) de soporte local que recibe el nombre de función de forma del campo continuo de 
desplazamientos u;(x,y,z) o velocidades ú;(x, y, 2), 


= k(x, y, z)Ük (11.5) 


Qe 


uj(x,y,z) = k(x, y,Z)Uk = ù;(x,y,z) 


Qe 


Qe Qe 
Esta función N;z(x,y,z), que actúa sobre un dominio acotado Qf denominado elemento 
finito, permite aproximar dentro de dicho dominio los campos de desplazamientos uz(x, y, z), 
velocidades úx(x,y,z) y aceleraciones ¡iz(x,y,z) mediante la valoración de sus respectivas 
magnitudes U;, Ux y Üx en un número finito de puntos, denominados nodos, pertenecientes al 
dominio del elemento finito (2%. De esta forma puede establecerse los campos derivados del 
desplazamiento, como lo es entre otros la deformación de Almansi ej. = vo up. Esto es, 


_wW,. 
= V? uj 
Qe 


= V? N; Uk 
Qe 


uj(x,y,z) = k(x, y, 2JUg (11.6) 


Qe 


= ei; 


Qe 


Qe 
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Se denomina método de los elementos finitos al procedimiento numérico que surge de 
utilizar esta aproximación polinómica para las funciones de campo [37]. Esta aproximación 
reduce las infinitas incógnitas de la función de campo a un número finito de incógnitas, defi- 
nidas en ciertos puntos preestablecidos como nodos del elemento finito. 


Sustituyendo la aproximación 11.5 y 11.6 en la ecuación 11.4, puede escribirse la ecuación 
de equilibrio de potencias a partir de la siguiente aproximación: 


J oV? Naav Us 
ye Qe 


= F tNudS+ | pbiNxdV — | PNuNjÜjdv Ù; 
S ye ye Qe 


Qe Qe 
(11.7) 

Pero esta ecuación se cumple para cualquier velocidad Ueliós por lo tanto la igualdad 
establecida en la ecuación 11.7 es independiente de esta velocidad, obteniéndose de aquí la 


siguiente ecuación de equilibrio de fuerzas dinámicas para el sólido discreto: 


Bijk oe MJ 

cr 
Jo V? N; dV = puNads+ | pbiNadV - f pNuNyav| mi pras), 
ye a ye ee. ye s e 
E > pi a KA 

pe öë J öt 


(11.8) 
3 int $ 
MRUj+ SE =f E QG 

siendo f?” | pe |] ayl | oe los conjuntos ordenados, en forma de matrices columna, de 
las fuerzas interna, másica y externa que se desarrollan en cada punto del sistema discreto 
que aproxima el continuo, Ü j | e la aceleración en dichos puntos, My; ğe la masa elemental, 
My; la matriz de masa y B;jk| oe = V¿N kl oe el tensor de compatibilidad de deformaciones o 
gradiente simétrico de la función de forma [23]. 


Desde un punto de vista mecánico-numérico, la no linealidad en la ecuación 11.8 puede 
estar originada por distintos fenómenos: 


= No linealidad constitutiva, que resulta de la pérdida de linealidad entre el campo de 
tensiones y deformaciones O;;— €;;, tal como ocurre en la plasticidad, daño, etc. Esta 
no linealidad ocurre debido al cambio de propiedades que sufre el material durante su 
comportamiento mecánico y se refleja en su tensor constitutivo C; jg. 


= No linealidad por grandes deformaciones, que es debida a la influencia no lineal 
que tiene el cambio de configuración del sólido en el campo de deformaciones. Este 
cambio de configuración también altera el tensor constitutivo C;jx1, y por ello estable- 
ce una relación no lineal entre tensiones y deformaciones. Además, estos cambios de 
configuración son producidos por grandes movimientos, traslaciones y rotaciones, que 
también producen cambios en el sistema de referencia local en los puntos del sólido, 
afectando por ello al tensor de compatibilidad de deformaciones B; jk 


a No linealidad por grandes desplazamientos, que a diferencia de las grandes deforma- 
ciones, sólo afecta al tensor de compatibilidad de deformaciones B; jk, porque en este 
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caso sólo ocurren cambios en el sistema de referencia local de los puntos del sólido 
como consecuencia de grandes movimientos. 


Estas posibles no linealidades pueden ocurrir todas a la vez o por separado. 


La ecuación 11.8 representa el equilibrio en el dominio un único elemento finito (2*, y 
su participación en dominio global Q se realiza a través del concepto de "ensamblaje" de 
esta ecuación de equilibrio junto a otras similares pertenecientes a otros elementos finitos 
que comparten el dominio del sólido. Esta operación se realiza mediante el operador lineal A 
que representa la suma entre las componentes de la fuerza, según corresponda a la posición y 
dirección de las contribuciones locales. 


En el caso que haya linealidad en el comportamiento del sólido, se cumple la siguiente 
relación de equilibrio global, cuya expresión resulta del ensamblaje de las ecuaciones de 
equilibrio local representadas en la ecuación 11.8: 


0=A [qna po 10 ne Afi! (11.9) 


La no linealidad en el comportamiento global del sólido se manifiesta como una fuerza 
residual A fal , provocada por el desequilibrio entre las fuerzas interiores J” | o» las fuerzas 
másicas f" | y las exteriores f” la Este desequilibrio, en un cierto instante de tiempo 
t del proceso dinámico, puede eliminarse mediante la linealización de esta fuerza residual 
A Ala (ecuación 11.9), en la vecindad del estado de equilibrio actual (i+ 1). Para ello es 
necesario forzar el equilibrio en el estado actual (i+ 1) y expresar dicha condición mediante 
una expansión en serie de Taylor truncada en su primera variación, 


; i i t TOAFA i t 
o= Aaso SA afata | AUS 


Qe 
. Ñ i ðÜ; ð Jor ð fint ðÜ; ð ea t f 
_ ihl t mi t , Ji k k J k „i41 t 
0="lAfdo = Afiat g moso - ðU, ` 9U;0U, 00, L [AU] 


(11.10) 


donde la aceleración y la velocidad deben expresarse mediante una aproximación lineal en 
diferencias finitas [23], como el método de Newmark. Sustituyendo en esta ecuación las 
fuerzas internas y másicas expresadas en la ecuación 11.8, se tiene, 


t 


E 7 i dU; 
S a 
a M,0,+ | oyy; Ni dV — cl TA J onun; av Pal 
ye Qe ye 
$. 
dOij des; S oO; o Ds OU, of 

1 I V?N dV LVN dV |= — 11.11 
ye des: ðU, dl t ye Ds Um dd ðU, ðU, f i 


-A HIJAU, Jh. 
A [AU Lo 


tal que particularizando esta ecuación de equilibrio dinámico para un material cuya ley cons- 
titutiva visco elasto-plástica es del tipo ©;j = p (DY (eij, p;)/de;;) = Cia: €; + Ena : Da 
para una relación cinemática del tipo e;; = Vu; = VN Uk y Dij = VEN ¡xÚx. resulta, 
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t 


0=4 M,0,+ | oy VENudV — $2" +A J onuni av 
ye Qe ye 


DU If 


+ f (VSN,,) CT, (VSN ye) av | + f (VEN) EZ (VSN ye) dV 


JU, ðU, | 
ye ye Qe 
-4 AU Joe (11.12) 
$ $ . T . 1 1 
0="(4f] + [Ad „i+ av, (11.13) 

Q Q Q 


donde Fis = El es el tensor de viscosidad tangente y J da = J" es el operador Jacobiano, que 
en un problema estático se convertiría en la matriz de rigidez tangente. Esta ecuación puede 
también presentarse en la siguiente forma matricial, donde se detallan los operadores que 


contribuyen a la definición del jacobiano: 


; tog ih oÙ Ù o], 1 
— il mi T T 141 
0= jar] = far| + Mp +K +05 E jao] (11.14) 
NMa 
To 


siendo esta última la ecuación de equilibrio linealizada, donde KĒ = A Sy. (VSN) TOE 
(VSN) dV representa la matriz de rigidez tangente, Mo = A Jye PN:NadV es la matriz de 
masa, D = A Sy. (VSN) ¡El : (VSN) dV es la matriz de amortiguación tangente, todas ellas 


definidas en todo el dominio Q; Ciy el tensor tangente correspondiente a la ley constitutiva 
utilizada en cada punto del sólido y f® = A | fse N : tdS+ fye PN : bdV] es la fuerza exte- 


. “q. 21" i t TE . . ., 
rior. La fuerza desequilibrada en el sólido *+! [Af] o Se elimina siguiendo una resolución por 
Newton-Raphson hasta que este residuo resulte despreciable, situación que se conoce como 
convergencia del proceso linealizado hacia la solución exacta (véase Figura 11.2). 


En la Figura 11.2 se describe el equilibrio espacial, dejando el tratamiento de la convergen- 
cia en el tiempo para ser abordado en un estudio más profundo de los métodos de resolución 
en el tiempo de la ecuación de equilibrio dinámico. 


11.4. Problemas estáticos y dinámicos no lineales 


A continuación se hace una breve presentación de los distintos tipos de comportamiento 
que introducen no linealidad e inducen al deterioro de una estructura. En principio, y en con- 
secuencia con el alcance de éste trabajo, se presentan brevemente los conceptos de elasticidad 
retardada y relajación de tensiones, que ocurren en problemas dinámicos y que dan lugar al 
denominado amortiguamiento viscoso, en los cuales la disipación de la energía depende de 
la velocidad. También la plasticidad y el daño, que contribuyen a la disipación de la energía 
independiente de la velocidad, situación que también ocurre en problemas estáticos. Aunque 
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Incremento de la carga en el instante: /=0+47 


Incremento: “n=” > 'lag, h =X I el: 
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Cálculo de la deformación 
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Gillae |T V¡N ¿Us Qe 
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Fuerza residual 


"IN, E A [s E BA” [am Tra e 


Verificación de equilibrio 
y convergencia 


¿"la ],>0? 


Nuevo incremento: 
HAr=> +1” 


Figura 11.2: Representación esquemática de la resolución del problema no lineal mediante Newton-Raphson. 
El campo de velocidades y aceleraciones deben calcularse según el método aproximación de la aceleración 
(véase a modo de ejemplo el método de Newmark) [23]. 


las grandes deformaciones también inducen no linealidad en el problema dinámico, no será 
tratada en este capítulo. 


Para determinar si un punto del sólido está en un estado de pequeñas deformaciones, se 
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debe cumplir que el determinante del gradiente de deformaciones sea cercano a la unidad: 


dV 


det |F| = dV, 


IR 


1 (11.15) 


Resultando en este caso la coincidencia entre las tensiones de Cauchy y de Piola Kirchoff 
Gij = Sij, entre la velocidad de deformación en la configuración actualizada y la correspon- 
diente magnitud infinitesimal D;; = €;;, y entre la densidad en las distintas configuraciones 
p = Po. Con estas condiciones particulares los desplazamientos y deformaciones son des- 
preciables frente a las dimensiones del sólido y por lo tanto puede escribirse la medida de 
deformación de la siguiente forma: 


ol oT os. 1|/0u duy 
€ => (FF -I)=V i= ORE (11.16) 


Además, para un material particular cuyo comportamiento es elástico y en pequeñas defor- 
maciones, ocurre la siguiente coincidencia en la definición de las deformaciones e = E = €, 
y la energía libre se escribe en la siguiente forma simplificada, 


1 
YW= —[le:C:eg€ 11.17 
I7 ) ( ) 


tal que sustituida en la segunda ley de la termodinámica (desigualdad de Claussius-Duhem) 
resulta la siguiente ley constitutiva [23], 


o=p03=C:e (11.18) 


donde el tensor constitutivo coincide exactamente con el obtenido mediante la ley de Hooke 
generalizada, y cuya expresión canónica es la siguiente: 


Cija = 28581 + U (Ex Ój1 + Suk) (11.19) 


siendo A y u las constantes de Lamé y 0;; es el tensor de Kroneker. El tensor de elasticidad 
de Hooke resulta definido positivo y posee las siguientes simetrías: 


Cijra = Cktij = Cijik = C jilk (11.20) 


Cauchy definía cuerpo elástico como "aquel en el cual las deformaciones en cualquier pun- 
to del sólido quedan determinadas por su estado de tensión y temperatura". En contraste con 
esta definición, se tendrá un material con comportamiento inelástico, por lo que es necesario 
establecer unas definiciones adicionales a las propias de la teoría de la elasticidad clásica, 
cuya formulación está relacionada con la historia del comportamiento del material. Esta si- 
tuación hace que no pueda garantizarse una relación biunívoca entre el tensor de tensiones y 
el de deformaciones, o dicho de otra forma, que no son relaciones invertibles una de otra. 


La influencia del tiempo produce en algunos sólidos comportamientos irrecuperables. Bá- 
sicamente pueden establecerse tres tipos de comportamientos no lineales en los materiales 
dependientes del tiempo: 


a Elasticidad retardada o "creep" donde ocurre crecimiento de deformación a tensión 


aplicada constante (véase Figura 11.3). 


= Relajación de tensiones, donde se produce pérdida de tensión mientras el nivel de 
deformaciones se mantiene constante. Este comportamiento, aunque no invertible, re- 


Problemas estáticos y dinámicos no lineales 


presenta la forma implícita inversa de la elasticidad retardada (Figura 11.3). 


a Viscoplasticidad cuyo comportamiento no lineal se debe a un crecimiento del cam- 
po de deformaciones inelásticas, pero esto ocurre siempre que el campo de tensiones 


supere unos umbrales preestablecidos (véase Figura 11.4). 


Problema de Creep o 
Elasticidad retardada 


- Tiempo de retardo t 


O0 


Modelo de Kelvin 


o'(t)= Ce (t) 


y” (0) — E È vis (0) 


o(t)=0"(t) +0" (t) = 
=C (0+Eé"(t) 
elt) = € (t) = e" (t) 


Problema de Relajación 


€0 


Modelo de Maxwell 


a(t) =C e(t) =0"(t) =5*(t) 


(= EE” 


Figura 11.3: Formas simplificadas de entender el comportamiento viscoso de Kelvin y Maxwell. 
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Problema de Viscoplasticidad 


Dispositivo de fricción 
lim 


[0] 
o” = Ce” 


lim 


¡o 
m Y 


Figura 11.4: Formas simplificadas de entender el comportamiento viscoplástico. 


Hay también materiales que tienen también mecanismos que conducen a comportamientos 
no lineales independientes del tiempo, que puede ser consecuencia de: 


= Plasticidad o comportamiento con flujo instantáneo. Este comportamiento puede 
matemáticamente establecerse como un caso particular del comportamiento viscoplás- 
tico, pero la física del problema es cualitativamente diferente (véase Figura 11.5). 


= Daño o degradación de rigidez, que produce en los materiales una pérdida de resisten- 
cia como consecuencia de una degradación en la elasticidad del material (Figura 11.5). 


Estos comportamientos pueden presentarse en forma aislada o participar todos ellos en dis- 
tinto grado. Sobre la modelización constitutiva y su influencia en el comportamiento estruc- 
turas se recomiendo profundizar el estudio en las fuentes orientadas a esta finalidad [11, 23]. 


11.5. Daño estructural. Caso sísmico teniendo en cuenta las 
curvas de capacidad y la demanda estructural 


Habiendo presentado en los apartados anteriores las bases para el estudio del daño mecá- 
nico, a continuación se presenta una formulación que permite obtener de forma sencilla una 
cuantificación sobre el daño estructural o deterioro máximo que puede alcanzar una estruc- 
tura durante el evento dinámico-sísmico sugerido por una normativa para una determinada 
zona. 


El análisis que se presenta a continuación se basa en suponer que el comportamiento es- 
tructura no lineal sigue los fundamentos de la teoría del daño mecánico [12, 22, 23, 25, 27, 


Daño estructural. Caso sísmico teniendo en cuenta las curvas de capacidad y la demanda estructural 


Problema de Plasticidad Problema de Daño 


aV 


o =(l-dCe=Ce -dCe 


Parte degradada 
del muelle 
o 0] 


O =C =C(e - e”) <= MAA --—> 
Dispositivo de fricción (1-d]C, ac 
C gim e 3} 


Figura 11.5: Formas simplificadas de entender el comportamiento elastoplástico y daño. 


29, 34, 36]. Esta teoría está basada en los principios mecánicos presentado en los apartados 
anteriores y cumple con los principios fundamentales de la termodinámica. 


No todos los materiales utilizados con fines estructurales siguen un comportamiento asi- 
milable al daño (degradación o pérdida de rigidez), pues hay algunos que se ajustan más a 
la plasticidad (desarrollo de deformaciones irrecuperables) y otro tienen una composición 
de comportamiento entre daño y plasticidad (degradación con deformaciones irrecuperables) 
(véase Figura 11.5). Dada una respuesta estructural, sólo se puede saber si ha tenido un com- 
portamiento de daño o plasticidad cuando se inicia la descarga, pues es el momento en que 
se diferencia la forma de las curvas de respuesta, mostrando en el caso de daño una descarga 
al origen con pérdida de rigidez, y en la plasticidad conservando la rigidez pero acumulando 
deformaciones irrecuperables. Mientras están en carga, la forma de la curva puede ser idénti- 
ca para daño y plasticidad. El hormigón armado resulta de un comportamiento mixto (daño y 
plasticidad), pero con preponderancia del efecto de degradación [12]. Esta afirmación puede 
ser confirmada al observar ensayos de laboratorio y puede ser cuantificada en base a la teoría 
de mezclas de sustancias simples [17, 19-21]. 


11.5 
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11.5.1. Daño global en función de la rigidez estructural 


La formulación que a continuación se presenta ha sido pensada con el objetivo de obtener 
información sobre el deterioro de estructuras, sometido a acciones estáticas y dinámicas (sís- 
mica) con muy poca información, cuyo origen hacen del procedimiento un método sencillo 
para una utilización rápida y eficaz. 


Suponiendo un comportamiento elástico de la estructura se hace un análisis tipo push-over 
y de aquí se obtiene la rigidez estructural inicial Kọ que se muestra en la Figura 11.6. 


Vv F , : 7 
Vz — 
yep 
| 
Respuesta ideal 
Elasto-plástica 
B 
Va — — — — 
Punto alcanzado / 
A utilizando la SINN ANS SN 
” máxima ductilidad DAAA 
Dr — y 
A K, > -A Rigidez con daño en vigas y pilares 
A, AKA A 


Figura 11.6: Cálculo de la rigidez estructural bajo cargas de push-over. 


Suponiendo conocido el máximo cortante en la base V que desarrollaría la estructura en 
el momento que se inicia su plastificación V, y suponiendo que el cálculo estructural se ha 
realizado una ductilidad u (obtenida de las normativas), conforme al tipo de la estructura 
en análisis, se tiene la siguiente expresión para el daño estructural máximo esperable (en el 
punto c, cuando desarrolla toda su ductilidad esperada), evaluado según la mecánica de daño 
continuo: 


Vy Y 
K, A, A 1 =1 
pl E E AT A AR (11.21) 
Ko Vy Y poo op 
A, Ay 


De aquí se desprende que el daño máximo, cuando la estructura desarrolla toda su ductilidad, 
es sólo función de la ductilidad adoptada. Así pues, podría decirse lo siguiente respecto del 
daño estructural máximo esperable, 


a Estructura dúctil: u=4 = De=0,75 


= Estructura frágil: u=2 = D.=0,50 
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que para llegar al comportamiento último de la estructura puede alcanzarse un nivel de daño 
mayor en una estructura dúctil que en una estructura frágil. Sin embargo, esta forma de medir 
el daño puede llevar a equívocos, pues podría también interpretarse lo contrario, es decir, que 
una estructura dúctil se daña un 25 % más que una frágil en su estado último. Esto obliga a 
objetivizar el cálculo de este índice de daño por comparación con el máximo daño que puede 
alcanzar la estructura. Así, el daño objetivo el < 1 alcanzado por una estructura en un punto 
p, se define como: 


DPS 2 = = y com OLDI (11.22) 


Se define ahora a p como el "Punto de Desempeño Estructural", resultante de la intersec- 
ción entre la curva de "Demanda Inelástica" y la curva de "Capacidad Estructural" (obtenida 
mediante un ensayo push-over). En estas condiciones, la ecuación 11.22 dará el máximo daño 
estructural que alcanzaría la estructura sometida al terremoto indicado por la normativa del 
lugar. Una comparación de este daño objetivo puede verse en las Figuras 11.7, 11.8 y 11.9 
que se muestran a continuación. 


1.8 , sl 
an /— Demanda elástica i 
[et] 
“a l5 a 
$ — 
1.2 A 4 
„Punto de comportamiento 
o9 L Ka Curva de capacidad 
0.6 H : Es 
L ~ Demanda inelástica 
0.3 ' 
0 | | | | | e | | | | | | | | | | | | | 
O 30 60 90 120 150 180 210 240 270 300 330 360 390 420 450 480 510 540 570 600 


S¿ (mm) 


Figura 11.7: Punto de Desempeño Estructural para forjados reticulares, sometido a un 
espectro inelástico de proyecto para suelos Tipo Il y u =2. 


11.5.2. Daño estructural objetivo 


Cuando de define una curva de repuesta ideal aproximada por tres tramos, en vez de los 
dos que se suele utilizar en el índice de Park [23], se obtiene una mejora sustancial en el 
cálculo del índice de daño estructural objetivo D%J, y por ello una mayor aproximación a 
ensayos experimentales y numéricos resueltos con elementos finitos. Así, la forma de las 
curvas quedará según se muestra en la Figura 11.10. 
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S, (8) 


Demanda elástica 
Punto de comportamiento 
C Curva de capacidad 
Demanda inelástica 


O 30 60 90 120 150 180 210 240 270 300 330 360 390 420 450 480 510 540 570 600 


S¿ (mm) 


Figura 11.8: Punto de Desempeño Estructural para edificios porticados con vigas planas, 
sometido a un espectro inelástico de proyecto para suelos Tipo II y u =2. 


Demanda elástica 


S, (g) 


Punto de comportamiento 
Curva de capacidad 


= 
— 


Demanda inelástica 


O 30 60 90 120 150 180 210 240 270 300 330 360 390 420 450 480 510 540 570 600 


S¿ (mm) 


Figura 11.9: Punto de Desempeño Estructural para edificios porticados con vigas de canto, 
sometido a un espectro inelástico de proyecto para suelos Tipo II y u = 4. 


Utilizando la nueva curva tri-lineal y la ecuación 11.22 se obtiene un daño estructural más 
ajustado al real para tres estructuras diferentes, y puede ser comparada con el análisis sim- 
plificado de Park y la técnica de elementos finitos, considerando la modelización constitutiva 
mencionada en la sección 11.2. Esta comparación puede verse en la Tabla 11.1. 
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INS 
Y 


Respuesta ideal 
Elasto-plástica 


Respuesta 


| 


A; | 


A 
E Estática no-lineal 


11.5 


Figura 11.10: Respuesta de una estructura a un análisis push-over. Respuesta no lineal real, 
respuesta ideal de Park y nueva representación tri-lineal de la respuesta. 


Estructura u =A4^u/Ay Ru= v,/ Vp a nn MEF pe 
Pórticos resistentes 5,15 1,88 E i a7 de o 
Pórticos vigas planas 1,55 1,34 a A P a a 
Forjados reticulados 1,57 1,94 E de 6 a a a 


Tabla 11.1: Comparación entre el daño estructural obtenido en forma simplificada utilizando la curva de Park, 
la nueva curva tri-lineal y el método de los elementos finitos (MEF). 


El análisis estructural previamente presentado permite la valoración objetiva del daño es- 
tructural en forma sencilla. Concretamente, la utilización de la ecuación 11.22 permite apro- 
ximarnos a una valoración muy cercana a aquellas que resultan de procedimientos de cálculo 
más costosos. Así pues, se puede conocer el nivel de daño estructural para la intersección de 
la curva de demanda con la de capacidad de la estructura. Para ello es necesario conocer: 


a La rigidez inicial Ko de la estructura, obtenida mediante un análisis elástico. 


= Obtención de una curva tri-lineal de capacidad estructural (Figura 11.10), a partir de la 
siguiente secuencia de puntos: 


e El punto A en la curva de capacidad, como el máximo cortante antes que se ini- 
cie el proceso no lineal en las vigas (cálculo de una estructura elástica - véase 
Figura 11.11). 
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e El punto B como el máximo cortante antes que se inicie el proceso no lineal en 
las columnas (cálculo de una estructura con pilares elásticos y rótulas en las vigas 
- véase Figura 11.11). 


e El punto C a partir de la ductilidad u recomendada por las normativas para el tipo 
de estructura que se está utilizando. 


= El punto P como intersección entre la curva de demanda y la curva de capacidad, de 
donde resultará Kp. 


= Cálculo del daño objetivo mediante la utilización de la ecuación 11.22. 


a) Aparición de rótulas plásticas en los b) Aparición de rótulas plásticas en los 
extremos de las vigas extremos de los pilares del primer nivel 


Figura 11.11: Estructura elástica hasta que aparecen las primeras rotulas plásticas en vigas (izquierda). 
Estructura con rótulas en todos los extremos de vigas (derecha). 


11.6. Simulación numérica de estructuras reforzadas con 
FRP y valoración del daño 


El uso de los materiales compuestos como materiales de refuerzo para estructuras ya cons- 
truidas es una de las tecnologías que están teniendo gran aplicación en la industria de la 
construcción debido a las ventajas que presentan [3]. Entre estas ventajas están las altas rela- 
ciones resistencia-peso y rigidez-peso, las cuales mejoran el comportamiento de la estructura 
existente sin alterar su configuración geométrica. Asimismo, son materiales livianos que no 
demandan cambios en la distribución del sistema estructural o en la cimentación. Además, 
para el caso de construcciones sometidas a ambientes especiales, presentan alta resistencia a 
la corrosión. 


No obstante, para optimizar el diseño estructural de los materiales compuestos en las es- 
tructuras, es necesario el análisis del comportamiento de los materiales compuestos donde 
se identifique la forma de participación de los diferentes componentes del compuesto y su 
interacción con otros materiales como el hormigón armado. De igual forma, se requiere el 
análisis del comportamiento global de la estructura, en el que se establezca la incidencia de 
los materiales compuestos como parte integral de los elementos estructurales reforzados. Un 
procedimiento eficiente para realizar estos análisis es la simulación numérica con elementos 
finitos. 


El comportamiento de las estructuras de hormigón armado reforzadas con FRP depende 
de múltiples parámetros, y el grado de influencia de muchos de ellos no ha sido aún deter- 
minado, lo que dificulta que experimentalmente se puedan estudiar todas las variables. Por 
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otra parte, la heterogeneidad y la anisotropía propias de los compuestos han hecho que no 
haya métodos de análisis sencillos y efectivos que permitan determinar su influencia sobre el 
comportamiento en las estructuras [19, 21]. 


Como consecuencia, el desarrollo para optimizar el uso del FRP en las estructuras se debe 
realizar de forma conjunta entre el campo experimental y el campo de la simulación numérica. 
No obstante, es necesario el desarrollo de nuevas herramientas y mejorar las existentes, en 
busca de que la simulación numérica a nivel global y local represente de una forma más 
aproximada el comportamiento real de las estructuras reforzadas con FRP, con la finalidad 
de obtener información que no se puede medir en los ensayos. 


Teniendo en cuenta que la simulación numérica en el análisis de los materiales compues- 
tos es un campo que está en desarrollo, y que su uso en obras civiles es relativamente nuevo, 
en este apartado se hace una presentación sintética de la teoría de mezclas serie/paralelo, que 
junto a los conceptos descritos en apartados anteriores, se convierte en una potente herra- 
mienta para el análisis numérico. Una descripción detallada de esta teoría se puede encontrar 
en [33]. 


11.6.1. Simulación del comportamiento de los materiales compuestos 


Los materiales compuestos están conformados por diferentes tipos de sustancias inorgá- 
nicas u orgánicas, cada uno de los materiales componentes teniendo su ley constitutiva que 
condiciona el comportamiento del conjunto en función de la proporción volumétrica y de la 
distribución morfológica que tenga dentro del compuesto. 


En la actualidad, gran parte de las investigaciones acerca del comportamiento de los mate- 
riales compuestos corresponden al campo experimental. Su estudio por medio de modelacio- 
nes numéricas está en desarrollo pero presenta algunas restricciones; incluso en simulaciones 
recientes de elementos reforzados con FRP, se representa el comportamiento constitutivo del 
material compuesto como el de un material homogéneo, sin tener en cuenta el aporte de sus 
componentes. 


Con el objetivo de solucionar esta limitación, se han propuesto diferentes teorías de si- 
mulación del comportamiento de los materiales compuestos como gestores del uso de los 
modelos constitutivos que representan el comportamiento de los materiales simples que com- 
ponen los materiales compuestos. Estas teorías de simulación son herramientas que pueden 
ser utilizadas dentro de una técnica de aproximación de elementos finitos para simular apro- 
piadamente el comportamiento de los materiales compuestos, entre ellas las más relevantes 
son [5]: 


= Teoría de la homogenización. Dentro de la mecánica de medios continuos se requiere 
del análisis desde dos escalas diferentes: una macroscópica en la que se determina el 
comportamiento de la estructura, y una microscópica en la que se obtiene el comporta- 
miento del compuesto partiendo de la respuesta de sus componentes. 


= Teoría de las mezclas. Considera que el comportamiento de cada componente define 
el comportamiento global del compuesto. Parte de la mecánica del continuo local bajo 
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el principio de interacción de las sustancias que componen el material, suponiendo que 
en el volumen infinitesimal del compuesto participan en conjunto todos sus compo- 
nentes. Asimismo, considera que cada uno de ellos contribuye al comportamiento del 
compuesto en la misma proporción que su participación volumétrica. 


Este apartado está orientado a la utilización de la teoría de mezclas, tomando como base 
la teoría de mezclas serie/paralelo propuesta por Rastellini [33], y que ha sido validada a 
través de la comparación de los resultados con diversas pruebas experimentales. Esta teoría 
se fundamenta en la teoría de mezclas clásica inicialmente estudiada por Trusdell y Topin 
[24]. Para más detalle sobre la simulación del comportamiento de los materiales compuestos, 
se recomienda recurrir a las siguientes fuentes [5, 17, 28]. 


11.6.2. Simulación numérica de un pórtico de hormigón armado refor- 
zado con CFRP 


El estudio de los efectos que han producido los grandes sismos en estructuras tipo pórtico 
pone en evidencia que las zonas más susceptibles a daño son las uniones viga columna y las 
bases de las columnas. Por ello, con el propósito de garantizar la estabilidad de las estruc- 
turas durante un evento extremo, en muchos estudios se hace énfasis en la necesidad de la 
rehabilitación y reparación de las estructuras antiguas o de las estructuras construidas antes 
de las actuales normas de diseño, siendo una de las alternativas de refuerzo el uso de FRP. 


Utilizando la teoría de mezclas serie/paralelo [14-17, 32, 33] en un programa de elemen- 
tos finitos, se muestra en este apartado un estudio de estructuras porticadas reforzadas con 
CFRP (Carbon Fiber Reinforced Polymer). Se realizó un análisis no lineal estático incre- 
mental (push-over analysis) de diez estructuras planas con una misma geometría. Cinco de 
ellas son de hormigón simple y las otras son de hormigón armado, con distintas configura- 
ciones de refuerzo utilizando CFRP. 


Aunque el CFRP se emplea como refuerzo en estructuras de hormigón armado o de acero, 
se estudió también el comportamiento de las estructuras de hormigón simple reforzadas con 
el propósito de analizar únicamente la influencia del refuerzo de CFRP sobre las estructuras. 


11.6.2.1. Geometría y configuraciones del refuerzo 


En la Figura 11.12 se indica la geometría de los diez modelos junto con las armaduras 
de las viga y de las columnas para los pórticos de hormigón armado. En la Figura 11.13 se 
presentan las configuraciones del refuerzo con CFRP, y en la Tabla 11.2 se indica la nomen- 
clatura utilizada para identificar los modelos. 


El refuerzo en las columnas de estos modelos corresponde a dos capas orientadas a 0° y 
a 90”, para tener en cuenta que la eficiencia del encamisado de la columna depende de las 
direcciones en que se coloque la fibra. Asimismo, cuatro de los diez pórticos analizados tienen 
CFRP en la base de las columnas, dado que según los resultados experimentales [30, 31] su 
nivel de influencia es notorio en la capacidad y en la ductilidad en las columnas. 
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Figura 11.12: Geometría tipo de las estructuras porticadas. 


ER HE 


SASF - ASF SAF - AF 
s 
(a) Portico sin CERP (b) Portico con CFRP en la unión (Q) Portico con CERP en la unión 
y en la base de los pilares 
140 uo 


1 : 


SAFV - AFV SAFCV - AFCV 


05 


(d) Portico con CFRP en la unión (e) Portico con CFRP en la unión, extensión 
y extensión en la viga en la viga y en la base de los pilares 


Figura 11.13: Modelos de los pórticos con las diferentes configuraciones de refuerzo [19]. 
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Tipo de pórtico Hormigón sin Hormigón con 
P P armadura de acero armadura de acero 
Pórtico sin CFRP (Figura 11.13a) SASF ASF 
Pórtico con CFRP en la unión (Figura 11.13b) SAF AF 
Pórtico con CFRP en la unión y en la 
base de la columna (Figura 11.13c) Pare AEL 
Pórtico con CFRP en la unión y extensión 
en la viga (Figura 11.13d) SA ANN 
Pórtico con CFRP en la unión, extensión en la viga SAFCV AFCV 


y en la base de la columna (Figura 11.13e) 


Tabla 11.2: Nomenclatura de los modelos de los pórticos. 


Con el fin de tener un mejor comportamiento de la viga ante un desplazamiento lateral, 
en los modelos con refuerzo se encamisó la viga en las zonas cercanas a las uniones viga- 
columna. Se seleccionaron dos longitudes de refuerzo en la viga para hacer una comparación 
del comportamiento cuando el refuerzo se extiende desde el borde de la unión dos y cuatro 
veces la altura efectiva de la viga. 


La teoría serie/paralelo permite asignar el hormigón y la resina polimérica como material 
de las matrices de los compuestos, mientras que el acero y la fibra de carbono constituyen el 
refuerzo. En la Figura 11.14 se indican los materiales compuestos asignados a los elementos 
estructurales de los modelos, dependiendo de si tienen o no armadura. Los porcentajes vo- 
lumétricos en los que participan los componentes en cada material compuesto se señalan en 
Tabla 11.3, donde se especifica la dirección de las fibras del compuesto considerada como 
parte de la anisotropía del compuesto [4]. Las propiedades asignadas a los materiales simples 
se muestran en la Tabla 11.4. 


11.6.2.2. Descripción de los modelos constitutivos aplicados a los materiales 


Los modelos analizados por control de desplazamientos tienen una malla de elementos 
finitos rectangulares de 4 nodos. En lo que respecta al análisis del comportamiento de los ma- 
teriales simples, para determinar el daño en el hormigón y la matriz polimérica del compuesto 
se aplicó el modelo de daño de Kachanov [22, 23]. El acero se consideró como un material 
elastoplástico utilizando el algoritmo Euler [23], mientras que la fibra de carbono se analizó 
como un material elástico y lineal. En todos los casos se ha partido de la hipótesis que no hay 
daño inicial en el pórtico, por lo que los ejemplos corresponden al caso de rehabilitación con 
CFRP. Además se supuso una adherencia inicial perfecta entre el refuerzo y el hormigón. 


11.6.2.3. Análisis del comportamiento de los modelos de hormigón en masa y de hor- 
migón armado reforzado 
(a) Resistencia y rigidez 


En la Figura 11.15 se observa que los diferentes modelos tienen una rigidez inicial similar, 
su comportamiento se mantiene lineal aproximadamente hasta los 0,0015 m; luego, según la 


Simulación numérica de estructuras reforzadas con FRP y valoración del daño 


DETALLE DE LA VIGA DETALLE DEL PILAR DETALLE DE LA UNIÓN 
EN HORMIGON EN HORMIGON EN HORMIGON 


Sección 
sin CREP 
Secció 
con CRFP 
Sección Sección 
sin CRFP con CRFP 
DETALLE DE LA VIGA DETALLE DEL PILAR DETALLE DE LA UNIÓN 
EN HORMIGON ARMADO EN HORMIGON ARMADO EN HORMIGON ARMADO 
Sección 
sin CREP 
LE IS H Secció 
LO SS H Seccion = RENYA NYS EN rr 
A? RIRS s, AA SS LaS 
Ka A A con CRFP NOR E 
Sección Sección bi NX o 
sin CRFP con CRFP OXX <SZEN 


MATERIALES COMPUESTOS MATERIALES COMPUESTOS MATERIALES COMPUESTOS 
EN LA VIGA EN EL PILAR EN LA UNION VIGA -PILAR 

Material 1 |] Material 1 Material 4 
Material 5 Material 2 Material 7 
Material 6 Material 3 = Material 9 
Material 7 Material 7 |] Material 10 
Material 12 Material 3 Ea Material 11 
Material 13 Material 9 Material 14 
Material 15 Material 14 Material 15 


Figura 11.14: Configuración de los materiales compuestos en los pórticos [18]. 


configuración del refuerzo del sistema, la pérdida de rigidez varía de acuerdo con el desarrollo 
del daño y de plasticidad. La resistencia máxima de cada caso se indica en la Tabla 11.5. 


En lo que respecta al comportamiento de las estructuras de hormigón simple, la resistencia 
de la estructura aumenta casi un 40% cuando se coloca CFRP en la base de las columnas. 
Por el contrario, en los casos con CFRP en la unión y en la viga (SAF y SAFV), el aumento 
de la resistencia del sistema es pequeño, un 5% y un 6%, respectivamente. Asimismo, se 
observa que los modelos con refuerzo adicional en la viga (SAFV y SAFVC) muestran un 
comportamiento dúctil después de alcanzar su capacidad máxima y, en los dos casos conser- 
van una pérdida de resistencia menor al 24%. Se aprecia como el modelo SAFC tiene una 
tendencia de pérdida de resistencia similar al modelo SAF, donde el refuerzo en las columnas 
incrementa la resistencia del sistema pero no incide en su ductilidad. 


Por otra parte, en los pórticos de hormigón armado, al reforzarse sólo la unión (AF) hay 
un pequeño incremento en la resistencia del orden del 7%. Cuando el refuerzo se coloca en 
la base de las columnas (AFC) o más allá de la zona de confinamiento en las vigas (AFV), el 
aumento de la resistencia llega a un 27%. Al reforzar tanto la base de las columnas como la 
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Material Matriz de Acero Acero Cercos Cercos Matriz CFRP 0  CFRP 90° 


compuesto hormigón long. vert. horiz. vert. polimérica horiz. vert. 
1 100,0 - - - - - - - 
2 98,70 - 1,10 0,20 - - - - 
3 97,53 - 1,10 0,17 - 0,40 0,40 0,40 
4 97,71 - 1,09 - - 0,40 0,40 0,40 
5 98,70 - - - 0,10 0,40 0,40 0,40 
6 99,90 - - - 0,10 - - - 
T 98,80 - - - - 0,40 0,40 0,40 
8 87,13 - 12,67 0,20 - - - - 
9 86,04 - 12,59 0,17 - 0,40 0,40 0,40 
10 81,51 4,84 12,45 - - 0,40 0,40 0,40 
11 92,87 4,84 1,09 - - 0,40 0,40 0,40 
12 93,86 4,87 - - 0,07 0,40 0,40 0,40 
13 95,00 4,90 - - 0,10 - z E 
14 - - - - - 34,00 66,00 - 
15 - - - - - 34,00 - 66,00 


Tabla 11.3: Composición porcentual de componentes en los materiales compuestos de los pórticos [18]. 


Material Criterio umbral E, = E, = E y Oe O; Ge G; 
de daño [MPa] [MPa] [MPa] [kPam] [kPa-m] 
Hormigón Mohr-Coulomb 2,5-10* 0,20 30,0 3,0 50,0 5,0 
Acero Von Mises 2,1-109 0,00 270 270 2000 2000 
Matriz polimérica Mohr-Coulomb 1,2-10 0,20 87,5 29,2 36,0 3,0 
Fibra de carbono Von Mises 1,5-10% 0,00 2300 2300 2000 2000 


Tabla 11.4: Propiedades mecánicas de los materiales componentes de los compuestos en los pórticos. 


viga (AFCV), la resistencia se incrementa un 50%. A diferencia de los modelos con hormigón 
simple, en todos los casos con armadura hay ductilidad independientemente del refuerzo. 
En los pórticos de hormigón armado, el CFRP como refuerzo, aumenta la capacidad de la 
estructura, dado que al disminuir el nivel de daño en el hormigón, posterga la plastificación 
del acero. 


Es relevante notar que el refuerzo no aumenta la rigidez lateral de los pórticos, lo que coin- 
cide con la observación de Tastani y Pantazopoulou [35]. Con el refuerzo hay un incremento 
de la resistencia a cortante y a flexión, no obstante, el nivel de incremento de resistencia de- 
pende de la configuración del refuerzo que se seleccione y de la armadura de los pórticos. En 
lo que respecta a la ductilidad, cabe anotar que su aumento depende de la distribución de la 
armadura; si la armadura es insuficiente, el refuerzo con CFRP incrementa considerablemen- 
te la ductilidad del sistema; pero si la armadura aporta por si misma ductilidad al sistema, al 
reforzarlo el incremento en la ductilidad es imperceptible. 


Simulación numérica de estructuras reforzadas con FRP y valoración del daño 


K"" =31764kN mo 


= 270004N / m 


0 0.005 0.01 0.015 0.02 0.025 0.03 00 0.04 
A (m) 
+- SASF —4—SAF —4-—SAFC —80-—SAFV ——SAFCV ——ASF —4- AF —4—AFC 40- AFV -—6-AFCV 


Figura 11.15: P— u pórticos de hormigón simple y de hormigón con armadura de acero [18]. 
Hormigón sin armadura SASF SAF SAFC SAFV*  SAFCV 


Fuerza máxima (kN) 62,36 65,48 87,65 66,06 89,73 
APimax (m) 0,0030 0,0030 0,0048 0,0031 0,0048 


Hormigón con armadura ASF AF AFC AFV AFCV 
Fuerza máxima (kN) 183,62 196,54 233,98 234,49 274,96 
APmax (m) 0,0354 0,0318 0,0312 0,0233 0,0222 


Tabla 11.5: Fuerzas máximas y desplazamientos en la mitad del claro de la viga. 


(b) Daño o deterioro estructural 


La teoría de mezclas permite obtener resultados por cada componente, entre los cuales se 
puede conocer el estado tensional, deformacional y el daño o deterioro estructural. A conti- 
nuación se muestra la evolución de la variable interna de daño o el endurecimiento plástico 
en el compuesto y en los materiales componentes. En la Figura 11.16 se muestra la distri- 
bución de daño en los pórticos para el estado de carga última. Los valores se recogen en la 
Tabla 11.6. 


En las estructuras de hormigón simple el daño se localiza y exhibe una considerable pérdi- 
da de capacidad, mientras que en las estructuras de hormigón armado el daño se distribuye a 
lo largo de los elementos estructurales y la pérdida de capacidad de la estructura es pequeña. 


Dependiendo del daño que se produce en los diferentes casos, se aprecia que en los pilares 
y en las vigas es necesario que la longitud del refuerzo externo sea suficientemente larga para 
evitar el daño localizado en el hormigón y retrasar la plastificación del acero, previniéndose 
el fallo prematuro del sistema. Como ejemplo, se observa que en las estructuras de hormigón 
reforzado -al comparar los modelos AFV con AF, y AFCV con AFC (véase Figura 11.13). 
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SASF Dc= 0,998 ASE Dc = 0,826 


Pa=1.49 KN Pw=183.61 KN 


Au=0.0375 m Ayyz0.0373 m 


Pa=2.09 KN P196.36 KN 


Au=0.0375 m Au=0.0375 m 


SAFC D¿= 0,997 AFC D¿=0,796 


FAM 


Pyy=2.63 KN Pyy=233.67 KN 


Au=0.0369 m Ayi=0.0366 m 


SAFV* D¿= 0,949 AFV Dc = 0,808 
Pu=51.40 KN Pyy=219.87 KN 


Ayz0.0143 m Ayi=0.0375 m 


SAFCY"  De= 0,933 


Ay=0.0114 m Ay=0.0371 m 


Figura 11.16: Daño estructural producido por un desplazamiento u en el extremo izquierdo. 
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Conclusiones l 11.7 


los pórticos con mayor longitud de refuerzo en las vigas tienen menor pérdida de rigidez, 
aunque presentan mayor daño en la unión y tienen una resistencia lateral un 12% superior. 
Asimismo, al contrastar los casos AF con AFC o SAF con SAFC, se observa que al colo- 
car el refuerzo en la base de las columnas, la zona con daño en toda la sección transversal 
se traslada del apoyo de la columna al borde en el que termina el CFRP, con lo cual, el re- 
fuerzo aumenta la resistencia del sistema. Utilizando la ecuación 11.21 se puede obtener el 
nivel de daño global que ha alcanzado cada una de las estructuras porticadas en el estado 
último de desplazamiento 4,1, = 0,0036 m al que se han exigido todas las estructuras (véanse 
Figura 11.6 y Figura 11.15). 


De éstas pruebas se confirma que al reforzar las estructuras con CFRP se aumenta la 
ductilidad y se reduce la pérdida de resistencia del sistema y el daño máximo estructural. Sin 
embargo, el nivel de eficiencia del refuerzo de compuesto depende tanto de la configuración 
de refuerzo que se seleccione, como del tipo y estado de la estructura original. 


Sin CFRP 
Modelo de Pórtico SASF' ASF* 
Pair [KN] 1,490 183,610 
Ult [m] 0,036 0,036 
Ko [kN/m] 27000 29382 
Ko [KN/m] = Pag ju ult 41,389 5100,278 
De = 1 — (K¿/K0) 0,998 0,826 
CFRP Unión Viga - Columna CFRP Unión Viga - Columna y Base 
Modelo de Pórtico SAFİ AF! SAFC' AFC* 
Pair [KN] 2,090 196,360 2,630 233,670 
uui [m] 0,036 0,036 0,036 0,036 
Ko [kN/m] 28000 31764,7 28000 31764,7 
Ko [KN/m] = Puit /u ult 58,056 5454,444 73,056 6490,833 
De = 1 — (K: /Ko) 0,998 0,828 0,997 0,796 


CFRP Unión Viga extendida 
Columna y Base 


CFRP Unión Viga extendida 
Columna y Base extendidas 


Modelo de Pórtico SAFVÝ AFV! SAFCV!* AFCVi 
Pair [KN] 51,400 219,870 68,780 259,650 
uus [m] 0,036 0,036 0,036 0,036 
Ko [kN/m] 28000 31764,7 28500 31764,7 
Ke [kN/m] = Pirja ult 1427,778 6107,5 1910,556 7212,5 
D.=1-(Ke/Ko) 0,949 0,808 0,933 0,773 


Y Sin Acero; * Con Acero. 


Tabla 11.6: Daño global, según ecuación 11.20, alcanzado por cada estructura para el desplazamiento último 
Unir = 0,0036 m. Máximo daño (0,998) en el modelo SASF y mínimo daño (0,773) en el modelo AFCV. 


11.7. Conclusiones 


La formulación presentada en este trabajo permite analizar numéricamente el comporta- 


miento de materiales compuestos y de sus componentes dentro de una estructura. Por una 
parte permite que cada componente del compuesto se analice utilizando la ecuación cons- 
titutiva que sea mas conveniente para predecir su comportamiento (elasticidad, plasticidad, 
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daño, etc.). Por tanto, esta formulación conduce a una herramienta numérica apropiada para 
el análisis de estructuras compuestas con diversos materiales con diferentes tipos de no linea- 
lidades. Al mismo tiempo, los costos computacionales son reducidos al utilizar un mallado 
continuo único para el análisis de la estructura y sus materiales componentes. 


El refuerzo con FRP en los pórticos de hormigón simple y armado, es una técnica viable 
para incrementar la capacidad de la estructura. Además, siempre que se configure adecua- 
damente el refuerzo en los elementos, se mejora el control de daño y se logra aumentar la 
ductilidad global del sistema. Según los resultados obtenidos en este análisis, tanto en los 
modelos de hormigón simple como en los de hormigón armado, el aumento de la resisten- 
cia global del sistema estructural depende de la configuración de refuerzo que se coloque. 
Por ello, es necesario que, a través de la simulación numérica utilizando la teoría de mezclas 
serie/paralelo, se optimicen los análisis de las estructuras con materiales compuestos para 
mejorar el diseño del refuerzo de los edificios y de los puentes que necesitan rehabilitación o 
reparación. 


La simulación numérica es una herramienta útil y de gran aplicación en el diseño de so- 
luciones de refuerzo de estructuras de hormigón armado con materiales compuestos. La si- 
mulación numérica con elementos finitos junto con la teoría serie/paralelo y los modelos 
constitutivos de los materiales componentes que ésta controla, permite analizar una estruc- 
tura con diferentes alternativas de refuerzo, involucrando los modelos constitutivos de los 
componentes de los materiales compuestos. Esto permite además que se pueda seleccionar 
las configuraciones más eficientes y se determine la capacidad de resistencia del sistema es- 
tructural reforzado, identificando el nivel de daño que se puede alcanzar. Esta aplicación es 
posible a dos escalas: una local, en donde se estudian los elementos estructurales, en la que 
los materiales compuestos se utilizan como refuerzo a cortante, a flexión o a compresión; y 
una escala global, en la que se analiza el comportamiento del sistema estructural al reforzar 
los elementos que lo requieran de acuerdo con las normativas existentes. 


11.8. Líneas futuras 


La formulación que se propone en este capítulo es bastante general, no obstante es necesa- 
rio extenderla a problemas termomecánicos, incluido el efecto del fuego sobre las estructuras. 
Esto abriría un campo de estudio muy importante que cubriría uno de los puntos más débiles 
y vulnerable de los materiales compuesto de matriz polimérica reforzada con fibras. 


También se está trabajando en la simulación numérica del comportamiento de compuestos 
con fibras cortas, abriendo así la posibilidad de capturar el incremento de ductilidad estructu- 
ral que aporta este tipo de refuerzo. 


Actualmente se está utilizando tejidos para reforzar estructuras. Esto ofrece una excelente 
alternativa a los laminados multidireccionales. La necesidad de tratar este problema e incor- 
porarlo a la formulación previamente presentada, es una necesidad por el gran crecimiento y 
aplicabilidad que están teniendo estos materiales compuestos. 


Por último, la combinación de la teoría de mezclas y homogeneización abren un potencial 
de análisis muy grande que permitiría abordar complejas situaciones que suelen presentarse 
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en las estructuras de materiales compuestos, como así también en el refuerzo de las mismas. 
Actualmente se está trabajando en esta línea con buenos resultados. 
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